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ABSTRACT 

Gas giants orbiting interior to the ice line are thought to have been displaced from their forma- 
tion locations by processes that remain debated. Here we uncover several new metallicity trends, 
which together may indicate that two competing mechanisms deliver close-in giant planets: gentle 
disk migration, operating in environments with a range of metallicities, and violent planet-planet 
gravitational interactions, primarily triggered in metal-rich systems in which multiple giant planets 
can form. First, we show with 99.6% confidence that giant planets with semi-major axes between 0.1 
and 1 AU orbiting metal-poor stars ([Fe/H]<0) are confined to lower eccentricities than those orbiting 
metal-rich stars. Second, we show with 98.3% confidence that eccentric proto-hot Jupiters undergoing 
tidal circularization primarily orbit metal-rich stars. Finally, we show that only metal-rich stars host 
a pile-up of hot Jupiters, helping account for the lack of such a pile-up in the overall Kepler sample. 
Migration caused by stellar perturbers (e.g. stellar Kozai) is unlikely to account for the trends. These 
trends further motivate follow-up theoretical work addressing which hot Jupiter migration theories 
can also produce the observed population of eccentric giant planets between 0.1 and 1 AU. 
Subject headings: planetary systems 



1. INTRODUCTION 

Approximately 1% of stars host hot Jupiters, 
ousted from their birthplaces to short-period orbits 
(|Wright et al.l I2012T) via mechanisms that remain de- 
bated. Proposed the ories fall into two classes: sm ooth 
disk migration (e.g. iGoldreich fe Tremaind 11980) . and 
migration via gravitational perturbations, either by stars 
(e.g. stellar binary Kozai, IWu fe Murravl 120031) or sib- 
ling p lanets (including planetary Kozai, e.g . iNaoz et al.l 
1201 It scatte ring, e.g. iRasio fc Fordl Il996t and secular 
chaos , e.g. IWu fc Lithwickl 1201 ID . (See iDawson et ail 
120131 DM J 13 hereafter, for additional references.) We 
consider the latter class as also encompassing gravita- 
tional perturbations p roceeded by disk migration (e.g. 
IGuillochon eT^Wmm . 

Migration processes must not only produce hot 
Jupiters — heavily studied, extensively observed gas gi- 
ants orbiting within 0.1 AU of their host stars — but 
also populate the region from 0.1 to 1 AU, interior to the 
ice line (beyond which large, rocky cores can grow and 
accrete) but outside the reach of tidal damping forces 
exerted by the host star. We call this semi-major axis 
range the "Valley," b ecause it roughly corresponds to the 
"Period Valley" (e.g. I Jones et al.ll2003l ). the observed dip 
in the giant planet orbital period (P) distribution from 
roughly 10 < P < 100 days. The Valley houses gas gi- 
ants both on highly eccentric and nearly circular orbits. 
Gas disk migration is unlikely to excite large eccentric- 
ities (e.g. iDunhill et al.l f20T3h whereas dynamical inter- 
actions are unlikely to produce a substantial population 
of circular orbits. Therefore this eccentricity distribu- 
tion may point toward intermixing between two differ- 
ent migration mechanisms, one gentle and one violent. 
Another orbital feature — the bimodal distribution of 
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spin-orbit alignments among hot Jupiters — is some- 
times interpreted as evidence for two migration mecha- 
nisms (iFabrvckv fe Winnll2009t [Morton fc : Johnsorl20lH 
INaoz et alJ 12012ft . However, it ma y result from stellar 
torques on the proto-planetary disk ^Batygin 2012fl. grav- 



ity w aves that misalign the star's spin axis (Roge rs et al 
2012 ) , or two regimes for t idal realignment ()Winn et al 



2010t lAlbrecht et al.ll2012l ). Because tides are negligible 



in the Valley (except at the most extreme periastron, 
e.g. HD-80606-b, HD-17156-b), we can interpret trends 
more easily. Excited inclinations and eccentricities can- 
not have been erased by tidal damping. 

If two common mechanisms indeed deliver close-in gi- 
ant planets, physical properties of the proto-planetary 
envir onment may d e termi n e whi ch is triggered. A decade 
ago, iSantos et al.l (|2001l . I2004T I discovered that giant 
planets more commonly orbit metal-rich stars, support- 
ing the core accretion formation theory. Independent 
and follow -up studies confirmed thi s trend for giant plan- 
ets (e.g. iFischer fc Valentil 120051 iSozzetti et all 12009 
Johnson et all I2010L iSousa et alTl201lL iMortier et al 




but not small planets (Ribas & Miralda-Escudc 



Buchhav e et al.1 20121). Ne i ther t he ISantos et al 



( 2004) nor the IFischer fc Valentil (|2005f ) samples exhib 



ited correlations between stellar metallicity and plane- 
tary period or eccentricity, but now the radial-velocity 
(RV) sample has quadrupled. It is time to revisit the 
planct-metallicity correlation, but now to gain insight 
into the dynamical evolution of planetary systems fol- 
lowing planet formation. 

Another motivation is the puzzlingly low occ urrence 
rate of hot Jupiters in the Kepler vs. RV sampl e (jYoudinl 
[20TTI : IHoward et all I2012T : IWright et all 12012) . Kepler 
targets have systematically lower metallicities than RV 
targets. We will show that differences in the planetary 
period distribution — not just the overall occurrence rate 
— between metal-rich and metal-poor stars may account 
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for the discrepancy. 

We uncover new stellar metallicity trends in the ec- 
centricities of giant Valley planets (§2), eccentricities of 
giant planets tidally circularizing (§3), and giant planet 
period distribution (§4). These correlations point toward 
planet-planet interactions as one of two mechanisms for 
delivering close-in gas giants (§5). 

2. ECCENTRIC VALLEY PLANETS ORBIT METAL-RICH 

STARS 

Valley gas g iants are unlikely to have formed in situ 
(Rahkov 2006) and exhibit a range of eccentricities (e) 
(Figure 1). Here we consider RV-detected giant plan- 
ets with msin i > 0.1Afj up , (queried from the Exo- 
planet Orbit DatabaseQ [EOD] . iWright et al.|[20Tlh . We 
impose stellar cuts of 4500K< T c s,* < 6500K and 
logg* > 4 throughout the paper to limit our sample to 
main-sequence stars with with reliable effective tempera- 
tures and surface gr avities in the Kepler Input Catalogue 
(IBrown et alJ[2T)llh . (We use the Kepler sample in §5.) 

Under the two migration mechanisms hypothesis, Val- 
ley planets on nearly circular orbits moved in smoothly 
through the gaseous proto-planetary disk, whereas those 
on eccentric orbits were displaced through multi-body 
interactions. In Figure 1, we emphasize planets with 
large eccentricities by plotting 1 — e 2 . This quantity 
is related to the specific orbital angular momentum, 
h = W a(l — e 2 ), an important parameter for dynami- 
cal interactions. This scale also minimizes eccentricity 
bias. For example, as a result of noise and eccentricity 
bias, a planet truly on a circular orbit could have a mea- 
sured e ~ 0.1. However, on this scale, e = 0.1 would be 
nearly indistinguishable from e = 0. 

We divide the sample into planets orbiting metal- 
rich stars ([Fe/H]>0, blue circles) vs. metal-poor stars 
([Fe/H]<0, red squares). Only the metal- rich stars host 
Valley planets with large eccentricities. The eccentrici- 
ties of these 68 planets extend up to 0.93. In contrast, 
the 26 Valley planets orbiting metal-poor stars are con- 
fined to low eccentricities (e < 0.43). Overall, 26% of 
Valley planets orbiting metal-rich stars have eccentrici- 
ties exceeding that of the most eccentric one orbiting a 
metal-poor star. 

We assess the statistical significance of the low eccen- 
tricities of Valley planets orbiting metal-poor stars. We 
perform a Kolmogorov-Smirnov (K-S) test on the null 
hypothesis that the eccentricities of the metal-rich and 
metal-poor sample are drawn from the same distribution. 
We reject the null hypothesis with 99.4% confidence. Us- 
ing a test more sensitive to the tails of distributions, 
Anderson-Darling (A-D), we reject the null hypothesis 
with 99.8% confidence. The results are insensitive to the 
exact metallicity cut and statistically significant for any 
cut located between -0.1 and 0.1 dex. Finally, we perform 
a Monte Carlo simulation, randomly selecting 26 planets 
from the total sample of 94 Valley giants. In only 0.4% of 

2 Five planets fulfilling our selection criteria have eccentricities 
fixed at in the EOD fits. We perform Monte Carlo Markov 
Chain fits to the RVs of 14-And-b (e = 0.026±q HD-81688-b 
(e = 0.03lt° ° 2 °), and Xi-Aql-b (e = 0.26 ±0.04) usine lSato eTaLl 
(t2008lV s data; adopt Ijohnson et al.l j20TllV s e = 0.03(< 0.28) for 
HD-96063-b; and remove HD-104067-b because the RVs are un- 
available. 



10 7 trials is the maximum eccentricity of the 26 planets 
less than or equal to the observed e = 0.43. Therefore, 
with 99.6% confidence, we reject the hypothesis that the 
confinement to low eccentricities of the planets orbiting 
metal-poor stars results from chance. Although the ex- 
act statistical significance is somewhat sensitive to the 
definition of the Valley, which defines the sample size — 
e.g., if we cut off the Valley at 0.8 AU instead of 1 AU, 
our confidence is 98.3% — it is evident in Figure 1 that 
the trend occurs t hroughout the Valley. 

As suggested bv Uohansen et all (|2012f > in the context 
of the mutual inclinations of Kepler multi-planet sys- 
tems, one might expect a threshold metallicity to trigger 
instability. Decreasing planets' semi-major axes (a) via 
gravitational perturbations requires interactions between 
at least two (and probably more) closely-spaced giant 
planets. It may be that only metal-rich proto-planetary 
environments can form such systems In contrast, plan- 
ets on circular orbits would have arrived via disk migra- 
tion, which can occur regardless of metallicity. 

We note that beyond 1 AU, the metal-rich and 
metal-poor sample have similar eccentricity distribu- 
tions. Planets with a > 1 AU have not necessarily 
changed their semi-major axes: they may have formed 
where we observe them. These planets on eccentric or- 
bits near their formation location may have exchanged 
angular momentum with another planet or star without 
requiring the abundance of closely-packed giant planets 
necessary to drastically alter a. 

3. PROTO-HOT JUPITERS ORBIT METAL-RICH STARS 

We turn to planets experiencing significant tidal dissi- 
pation, detected^ by non-Kepler transit surveys (Figure 
2) and followed up with RV measuremen ts. We use the 
stellar and planetary cuts described in §1. ISocrates et al.l 
(|2012|) and DM J 13 used this sample to calculate the 
abundance of moderately-eccentric proto-hot Jupitcrs. 
Advantageously for this sample, transit surveys are less 
inclined to target metal-rich stars, yielding planets or- 
biting jnetaljTJO^r^arefor comparison. To be consistent 
with ISocrates et ail (|2012[ ) and DMJ13 and to avoid ec- 
centricity bias, we classify planets with e > 0.2 as eccen- 
tric. 

The striped region contains planets undergoing tidal 
circulariz ation along tracks of constant angular momen- 
tum (see ISocrates et all 120121 DMJ13) to final orbital 
periods Pfmai between 3 and 10 days. (The traditional 
boundary for hot Jupiters is 10 days, and 3 days is the 

3 RV systems containing multiple known giant planets do ap- 
pear to have systematically higher metallicities than those con- 
taining one, but the statistical significance is marginal. We note 
that planets may be scattered to distancesbeyond current RV de- 
tection or ejected, so systems with only one known giant planet 
perhaps originally had more. 

4 Some planets have e fixed at in EOD fits. We remove those 
with poorly-constrained eccentricities: CoRoT-7-b, HAT-P-9-b, 
OGLE-TR-10-b, OGLE-TR-lll-b, TrES-l-b, TrES-4-b, WASP- 
13-b, WASP-39-b, WASP-58-b, XO-l-b, XO-5-b. We include plan- 
ets whose eccentricities are constrained to be small (e < 0.2), by 
our fits (CoRoT-13-b, CoRoT-17-b, WASP-16-b) or the literature 
(CoRoT-7-b, HAT-P-l-b, HAT-P-4-b, HAT-P-8-b, HAT-P-12-b, 
HAT-P-27-b, HAT-P-39-b, OGLE-TR-211-b, KELT-2-Ab, WASP- 
7, WASP-ll-b, WASP-15-b, WASP-21-b, WASP-25-b, WASP-31- 
b, WASP-35-b, WASP-37-b, WASP-41-b, WASP-42-b, WASP-47- 
b, WASP-61-b, WASP-62-b, WASP-63-b, WASP-67-b). See the 
EOD for each planet's orbital reference. 
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Fig. 1. — Left: Valley (gray region) giant planets orbiting metal-rich stars ([Fe/H]>0, blue circles) have a range of eccentricities; 
those orbiting metal-poor stars ([Fe/H]<0, red squares) are confined to low eccentricities. For reference, above the dashed line (a tidal 
circularization track ending at 0.1 AU) planets are unlikely to experience significant tidal circularization. We plot the quantity 1 — e 2 to 
emphasize high-eccentricity planets. Right: Eccentricity distributions of Valley planets orbiting metal-rich (blue solid) and metal-poor (red 
dashed) stars. 




Fig. 2. — Left: Giant planets discovered by non-Kepler transit surveys, orbiting metal-rich (blue circles) and metal-poor (red squares) 
stars. The striped region encloses planets undergoing tidal circularization to 3 < Pfl na i < 10 days. Planets below the dotted line have 
e > 0.2, all of which orbit metal-rich stars. Right: Distribution of host star metallicities for planets in the striped region (left) with e >0.2 
(dotted line) and e < 0.2 (solid line). 



limit above which we still see eccentric giant planets. 
Those with -P nna i < 3 days have much faster tidal cir- 
cularization rates.) We only observe eccentric planets 
orbiting metal-rich stars (blue circles). We suggest that 
only giant planets forming in metal-rich systems with 
multiple giant planets are likely to be scattered onto ec- 
centric orbits that bring them clos e enough to the star to 
undergo tidal circularization (e.g. iFord fc Rasiol l2006). 

To assess the statistical significance, we perform a 
Monte Carlo simulation. We randomly select 8 host star 
metallicities from the 53 total host stars in the striped 
region. In only 2.3% of 10 7 trials was [Fe/H]> for 
all eight. We also perform a K-S (A-D) test, rejecting 



with 99.1% (99.2%) confidence the null hypothesis that 
the host star metallicities of planets in the striped re- 
gion with e >0.2 are drawn from the same distribution 
as those with e < 0.2. 

4. THE SHORT-PERIOD PILE-UP IS A FEATURE OF 
METAL-RICH STARS 

IHoward et all (|2012T > found a surprisingly low Kepler 
hot Jupiter occurrence rate (/hj, Kepler) — the expected 
number of giant planets per star with P < 10 days - 
compared to RV sur veys (/h.t.rv), a trend con firmed by 
iWright et all (|2012f ) and iFressin et al.l (|2013[ ); all sug- 
gested that the systematically lower metallicities of Ke- 
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Fig. 3. — Red striped: number of transiting giant planets de- 
tected by Kepler. Black dashed: expected number based on the 
RV-detected samplcEl The gray error bars are from uncertain- 
ties in Cnorm, not the Poisson uncertainties of each individual bin. 
The two distributions are consistent at long periods, but the Kepler 
sample lacks a short period pile-up. 

pier host stars may contribute to the discrepancy. In 
Figure 3, we compare the period distribution of tran- 
siting g iant planet candid ates dete cted by the Kepler 
survey (iBurke et al l 120131: see also iBorucki et ail 120111 
and IBatalha et al.l 120131 ) — applying a radius cut of 
8< -Rpianet < 20i? cart h — to that expected from the RV 
sample^ using a normalization constant C norm (defined 
below). The two distributions appear consistent beyond 
10 days but differ strikingly at short orbital periods: the 
Kepler period distribution lacks a short-period pile-up 
(in fact, the absolute Kepler giant planet occurr ence de- 
clines to ward short orbital per iods, as modeled bv lYoudir] 
120111 and iHoward et aTll2012l ) . 

Although Kepler Input Catalog ( KIC) metallicity es - 
timates are known to be uncertain (Bro wn et al.ll20lTI ). 
we can roughly divide the Kepler sample into metal-rich 
([Fe/H]>0) and metal-poor ([Fe/H]<0). In Figure 4, we 
compare the period distributions for Kepler giant planets 

5 The RV sample is not uniform; we plot it for qualitative 
comparison. The expect ed distribution d erived from the period 
distribution reported by lUumming et al. (2008) appears similar. 
We therefore interpret the short-period pile-up as real, not due 
to preferential detection. F or the quantitative calcu lations in this 
section, we use the uniform Fischer & Valcnti (2005) sample. 
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Fig. 4. — Number of transiting giant planets observed by Kepler 
without a stellar metallically cut (top), with [Fe/H]>0 (middle), 
and with [Fe/H]<0 (bottom). In the metal-rich sample (middle), 
we recover the shape of the short-period pile-up seen in the RV 
sample (black-dashed line, Figure 3). In contrast, the metal-poor 
sample (bottom) is depleted in short-period giants. 

orbiting metal-rich vs. metal-poor stars. When we limit 
the sample to [Fe/H]> (row 2), we recover the missing 
short-period pile-up, which the metal-poor sample (row 
3) lacks. Performing a K-S test, we reject with 99.94% 
confidence the hypothesis that the metal-rich sample and 
metal-poor sample are drawn from the same distribution. 
The results are insensitive to the exact metallicity cut. 

We compare the Kepler metal-rich(poor) sample to the 
RV metal-rich(poor) sample in Figure 5. In Figures 3 and 
5, we compare the observed number of transiting Kepler 
giant planets (red striped) to the number expected (black 
dashed) based on the RV sample, 



RV, trans 



= c n 



^Rvprob,. 



where -/Vrv is the observed number of RV planets per bin 
and prob trans (P) is the transit probability. We set the 
normalization constant, C n0 rm, using the values (com- 
puted below) of / H j,Ko P icr and / H j,rv: 



f todays at 

JHJ.RV 2^P=0 JV tre 



,Kep(P)/prob trans (P) 



/hj,k 



cplcr 



Ep=o yS ^Rv(P) 



Each error bar is due to the uncertainty in 
/h.t. rv//h.i. Kepler- 

To compute /Hj,Ko P ior, we fol- 
low IHoward et al.l (|2012l) . using our own stellar and 
planet ary cuts and the l atest sample of Kepler candi- 
dates (|Burke et al.l 120131 ). The Barbara A. Mikulski 
Archive for Space Telescopes (MAST) supplied the 
stellar parameters and the NExSci Exoplanet Archive 
the transit shape parameters (duration, depth, a/R+, 



RpUnct/R*)- We obtain^ / H j, 



Kepler 



= 0.37 



-0.07 



% 



6 We estimate the occurrence rates and uncertainties 
Markov Chain Monte Carlo using the Poisson likelihood. 
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Fig. 5. — Same as Figure 3 but for metal-rich (left) and metal-poor (right) subsamples. Left: Metal-rich Kepler sample (red striped) 
exhibits a short-period pile-up, but falls below RV expectations in the 3-5 day bin. Right: Metal-poor Kepler sample is not inconsistent 
with the metal-poor RV sample, but the latter is difficult to characterize due to small numbers. 

for giant planets with P <10 days (consistent with 
Howard et afll201~2l and iFressin et al.ll20Tl . 1.15±g;|l% 
for the metal-rich sample, and 0.25^0 07% for the metal- 
poor sample. To compute /hj,rv, we use the stellar 
and planetar y sample from the i conic planet-metallicity 



correlation ^Fischer fc Valentil I2005I) and associated 
stellar parameters ( Valenti fc Fischer! [2005), the last 
RV target list to be publicly released. We obtain 
/hj.rv = l-24±0 35% for giant planets with P <10 days 
(in agreement with iWright et alj|2012j ), 2.11±^|% for 



+0.39 
-0.19 



those orbiting stars with [Fe/H]>0, and 0.26 
for [Fe/H]<0. With no metallicity cut, /Hj.Kcpier is 
inconsistent with /hj,rv at the 2.4cr level. 

In the metal-rich comparison (Figure 5, left), we see 
greater consistency between the Kepler and RV distribu- 
tion than in the full sample (Figure 3). The metal-rich 
Kepler sample exhibits a short-period pile-up; the dis- 
crepancy between /Hj,Ke P ior vs. /hj,rv is now only lAa, 
with the greatest discrepancy in the 3-5 day bin. This 
improvement motivates a detailed follow-up analysis, in- 
cluding a more precise estimate of /hj,rv using the lat- 
est RV target lists. If follow-up studies find a significant 
discrepancy between the metal-rich Kepler and radial 
velocity samples, it could be due to the KIC metallicity 
estimates. Using spectroscopic metallicity measurements 



bv lBuchhave et all (|2012D . we find that high KIC metal- 
licities do correspond to high spectroscopic metallicities 
(with a scatter of about 0.2 dex), but the spectroscopic 
metallici ties are systematica lly higher by 0.1 dex, as also 
noted bv I Brown et al.l (|2011l ). However, we attribute the 
systematic offset to the fact that stars targeted for spec- 
troscopic follow-up are bright, main-sequence stars in our 
solar neighborhood and thus have systematically higher 
metallicities; in contrast, the KIC metallicities were com- 
puted assuming a low-metallicity prior, due to the Kepler 
targets being above the galactic plane. The planetary 
radius cut may also contribute to the discrepancy. The 
8-Rearth cut for the Kepler sample corresponds to the RV 
cut of mshu = 0.1Mj U p for a planet mad e of pure hy- 
droge n at a low effective temperature (e.g. iSeager et al.l 
|2007| ). However, close-in, low- mass planets may be in- 
flated to > 8i? oar th and may have a different period dis- 
tribution, contaminating the sample. In the metal-poor 
comparison (Figure 5, right), the Kepler and RV distri- 
butions do not appear inconsistent, but it is difficult to 
judge given the very small sample of RV-detected planets 
orbiting metal-poor stars. 

5. CONCLUSION 

We found three ways in which the properties of hot 
Jupiters and Valley giants depend on host star metallic- 
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ity: 



1. Gas giants with a < 1AU orbiting metal-rich stars 
have a range of eccentricities, whereas those orbit- 
ing metal-poor stars are restricted to lower eccen- 
tricities. 

2. Only metal-rich stars host eccentric proto-hot 
Jupiters undergoing tidal circularization. 

3. The pile-up of short-period giant planets, missing 
in the Kepler sample, is a feature of metal-rich 
stars and is largely recovered for giants orbiting 
metal-rich Kepler host stars. 

Hot Jupiters and Valley giants are both thought to 
have been displaced from their birthplaces. Therefore 
these metallicity trends can be understood if smooth disk 
migration and planet-planet scattering both contribute 
to the early evolution of systems of giant planets. We ex- 
pect disk migration could occur in any system, but only 
systems packed with giant planets - which most easily 
form around metal-rich stars - can scatter giant planets 
inward to large eccentricities (Trend 1), some of which 
tides shrink and circularize (Trend 2), creating a pile- up 
of short-period giants (Trend 3). Moreover, these trends 
support planet-planet interactions (e.g. scattering, secu- 
lar chaos, or Kozai) as the dynamical migration mecha- 
nism of close-in giant planets, rather than stellar Kozai. 
This is consistent with previous work by DMJ13 arguing 
that stellar Kozai does not produce most hot Jupiters, 
based on the lack of super-eccentric proto-hot Jupiters. 
We would not expect planet-planet scattering to typi- 
cally result in nearby companions to hot J upiters, which 
have b een ruled o ut in the Kepler sam ple bv lSteffen et al.l 
((20121) . fSee also lLatham et al.ll27mn 

One possible challenge for our interpretation is the 
lack of apparent correlation between spin-orbit misalign- 
ment and metallicity. However, spin-orbit misalignments 
are not necessarily caused by dynamical perturbations, 
and their interpretation is complicated because measure- 
ments have primarily been performed for close-in planets 
subject to tidal realignment. We recommend spin-orbit 



alignment measurements, via spectrosco py ((McLaughlin 
11924 iRossiterl 119241: [Queloz et all 120001) or pho tometry 
(jNutzman et al.ll201ll:lSanchis-Oieda et al.ll20Hl) . of Ke- 
pler candidates in the Valley, which are typically too 
distant to be tidally realigned. 

To support or rule-out the interpretation that these 
metallicity trends are signatures of planet-planet inter- 
actions, we further recommend: 1) theoretical assess- 
ments of whether planet-planet interaction mechanisms 
designed to account for hot Jupiters can simultaneously 
produce the observed population of eccentric Valley plan- 
ets, and 2) more sophisticated assessments of the trends 
we report here, using the target lists of recen t RV sur- 
veys and, as undertaken bv lFressin et al.l (12013D . a careful 
treatment of Kepler false-positives and detection thresh- 
olds. 
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